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Recently L.N. Bulaevskii, M. J. Graf and V. G. Kogan 
published two papers on the theoretical description of ex- 
periments on NbN thin film-based photon detectors Q ■ 
The central part of both works constituted the calcula- 
tion of the voltage produced by thermally-activated hop- 
ping of vortices across a thin film strip. The authors of 
Ref. 0, [I] used the approach developed in Ref. 0] , which, 
in turn, was based on the earlier work by Ambegaokar 
et al. |4[ describing the dynamics of a vortex in a film in 
terms of the Langevin equation rjx + U'(x) — C(t)- Here x 
is the position of the vortex across the strip (0 < x < w), 
((t) is the thermal noise source, and ij is the viscous drag 
coefficient. The energy U (x) of a vortex in a film of the 
width w < A = X 2 /d and the thickness d in the presence 
of the dc current / and the magnetic field H perpendic- 
ular to the film is given by [3| 



U(x) = eln 
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where e = 0q/167t 2 A, ^ » 0.34£, £ is the coherence 
length, A is the London penetration depth, (f> is the 
flux quantum, and c is the speed of light. The volt- 
age V(T, H, I) caused by uncorrelated jumps of vortices 
is calculated using the Fokker-Planck equation for the 
probability density f{x,t) 0, 



f = d x {U'(x)f + Tf] 



(2) 



The only difference between Eqs. (fT]) and $2$ used in 
Refs. 0, SI and Ref. is that the authors of Ref. 0, 
took ft = f/2 in Eq.©, while in Ref. this quantity 
was £i = 0.34£. The origin and the consequences of this 
discrepancy is discussed below. 

The results of Refs. 0,[|[ essentially reproduce those of 
Ref. I3j, however the vortex hopping rate R v derived in 
Ref. [3( differs by the factor F = 2 v (y — l)w/(, from that 
of Ref. 00]. The authors of Refs. 0, ij asserted that, 
for^ = e/T = 110 used in Ref. to fit the experimental 
data, R v of Ref. Q was overestimated by the factor 
F = T{v - \)wji ~ 3.5 x 10 36 because: 

1. The factor 2" in F comes from the factor 2 under 
the logarithm in Eq. (fT]) which, according to Ref. 0, Q, 
was missed in Ref. [3]. 

2. The factor v — 1 in F results from the use of the 



periodic boundary conditions for Eq. ([2]) in Ref. H as 
opposed to a 'realistic' boundary condition of Ref. [lj. 

3. The factor u>/£ in F results from the vortex interac- 
tion which, according to Ref. 0,[2[, leads to the statistical 
weight of a vortex P ~ L/w, as opposed to P ~ L/£ used 
in [3j, where L is the length of the strip. 

In this Comment we show that these statements are 
incorrect because they result from model artifacts of 
Refs. 0, Hj]. Below we address the issues taken into ac- 
count in Ref. [3] but neglected in Refs. 0, Q and dis- 
cuss their importance for a more consistent theory of 
thermally-activated hoping of vortices in thin films and 
the interpretation of experimental data. 



1. Core contribution 

The authors of Ref. apparently overlooked that 
£i = 0.34£ in Eq. JJ) of Ref. [1 absorbs both the 
factor 2 under the logarithm STa, 0] and the vortex 
core energy disregarded in Ref. 0, l2j . Here Eq . (fTJ re- 
sults from U(x) = e[ln[(2u;/7r£) s\n(-Kx/w)}+ f$\ [6j, where 
j3 = 0.38 accounts for the vortex core energy, so that 
^ = 0.34£ = e-^/2 in Eq.©, unlike & = £/2 used 
in 0,0. The value /? = 0.38 was extracted by comparing 
the lower critical field H cl = (0 o /47rA 2 )[ln(A/£) + 0.4971 
calculated from the Ginzburg-Landau (GL) theory [gj 
with the London result H cl = ((?!) /47rA 2 )[i^o(£/A) + 0\, 
where ft is the core contribution (associated with the spa- 
tial variation of the order parameter) not accounted for 
by the London cutoff, and K (^/X) » ln(2A/£) - 0.577 
for A 3> £. Matching these formulas for H c \ yields 
ft = 0.497 + 0.577 - In 2 = 0.38. In the London limit, 
the core line energy is independent of the sample geom- 
etry so ft w 0.38 is the same both in bulk samples and 
films with w 3> £, except for vortices spaced by x ~ £ 
from the surface. Taking /? into account significantly de- 
creases Rv(ft) — R v e~P' i l T as compared to R v calculated 
without the core contribution, while the variation of the 
core energy at the film edge affects the pre-exponential 
factor in Rv, as will be discussed below. 

One can see that the core contribution ft — 0.38 is 
no less essential than the London numerical correction 
ln(2/7r) = —0.45 in U(x), so taking ft into account is 
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important when comparing the model [lj, l2j with ex- 
periment. Indeed, neglecting the vortex core energy in 
Ref. [H overestimates R v by ~ exp(/3e/T) - 10 18 for 
e/T = 110. The importance of the vortex core contribu- 
tion in the vortex-related dynamic phenomena has been 
extensively discussed in the literature (see, for example, 
the recent work [qJ on the effect of the vortex core energy 
on the Berezinskii-Kosterlitz-Thouless transition). 

2. Boundary condition 

Here we show that the extra factor ~ v^ 1 in R v of 
Ref. [l[ does not come from the different boundary con- 
ditions used in Refs. l] and [3j, but rather from artifacts 
of the model of Ref. [1| . The vortex crossing rate was ob- 
tained in Ref. [l[ from the standard formula for a particle 
hopping between two potential wells |5| : 

R- 1 = D / e~ u ^' T dx / e u ^/ T dx, (3) 

JO Jxq 

where D = T/rj, xq ~ £, and X\ is a length smaller than 
x m at which U(x) is maximum. 

The authors of Ref. [l| assumed a model form of 
Ubgk(x) in Eq. Ubgk — U[x) where U{x) is given 
by Eq. © with ft = f/2 for x > x ~ ft U BG k{x) = 
for < x < xq, and Ubgk{x) = oo at x = 0. The in- 
finite repulsive barrier at the film edge was introduced 
artificially to trap vortices in the film by imposing the 
boundary condition of zero probability current S = if 
at x — for Eq. (0 . Vortex hopping in this model occurs 
as a 'pre vortex' [l[ is somehow placed in the film past 
this barrier, but it is unclear how this model can describe 
penetration of vortices in the film. 

The postulated form of Ubgk(x) — at < x < 
£ significantly overestimates the first integral Z = 
J Q Xl exp[—U(x)/T]dx in Eq. ©. To see how it happens, it 
is instructive to juxtapose Ubgk{x) with U(x) obtained 
by numerical simulations of vortices using the GL equa- 
tions, which take into account the vortex core energy and 
its change near the edge. These calculations have shown 
that the energy of a vortex, U (x) = (b + ax/^e, increases 
linearly with the distance x of the core phase singular- 
ity from the film edge up to x ~ £ [Hi [ll[ . This gives 
rise to a constant force ae/£ caused by a "string" of the 
suppressed order parameter between the core and the sur- 
face, where a ~ 0.1 — 0.3 and the constant b ~ 0.05 — 0.1 
accounts for the fact that U(x) > even at x — > due to 
local superconductivity suppression around a vortex core 
as it emerges from the film edge (lQ |. These features are 
essential for the evaluation of R v if U(x) > T at x < ft 

Substituting U(x) = (b + ax/t;)e in Z = 
J Q Xl exp[-[/ (x)/T]dx yields Z = ie~ hv /av for e~ al/ < 1. 
As follows from Eqs. ((T|) and ((3]), the factor e~ hv can be 
combined with eP v from the second integral in Eq. ([3]), 



so that the effect of the vortex core on the hopping rate 
R v ~ R v exp[(b — P)v\ is determined by the difference of 
core energies in the bulk and at the film edge. Here both 
a and b appear to be dependent of current [1 Of ] , indicating 
that the London notion of the rigid vortex core becomes 
hardly adequate at x ~ £. 

The calculation of Z ~ J^i^/xYdx ~ ^/v in Ref. Q 
is qualitatively consistent with the GL result. Here the 
cutoff ~ ^ where the London theory becomes invalid was 
used, and the upper limit can be extended to infinity 
if e~ av < 1 and I < I d , where I d = c(f> Q / '16tt 2 A£ is 
of the order of the depairing current. By contrast, the 
potential, U bgk{x) = at < x < xq, yields Z = xq ~ £ 
which overestimates Z by the factor ~ !/ > 1 as 
compared to both the GL results and Ref. [H . Treating a 
vortex like a particle in the London model combined with 
the Fokker-Plank equation does bring uncertain factors 
~ 1 in Z coming from the edge effects discussed above. 
Yet the simplified model of Ref. lj does not capture the 
qualitative behavior of Z ~ which follows from the 
more consistent GL theory and the approach of Ref. Q 
(also adopted in Ref. 0). A more realistic model of the 
vortex core penetration would require solving the time- 
dependent GL equations 

The above consideration shows that the claim of 
Ref. [l[ that the extra factor ~ v in i?" 1 comes from 
the 'realistic' boundary conditions as opposed to the pe- 
riodic U(x) of Ref. [3] is misleading. In fact, the solution 
of Eq. (O adopted in Ref. [3j is only defined inside the film 
< x < w and does not require any unphysical barriers 
at the film edges. Here the use of periodic U(x) in Eq. ([2]) 
is a standard method of satisfying the boundary condi- 
tions of a fixed probability flux S of vortices entering 
and exiting the film, which is equivalent to the method 
of images for solving the Laplace or diffusion equations. 
For example, Eq. (TTJ) can be obtained by either finding a 
proper analytical function or summing up potentials of 
an infinite chain of vortex-antivortex images outside the 
film. Moreover, if only the forward jumps of vortices are 
taken into account in the limit of e~ v C 1 1], Eq. © 
reduces to Eq. (7) of Ref. [|. This is not surprising 
because the exponentially small probability current S is 
mostly determined here by narrow vicinities of neighbor- 
ing minimum and maximum of U(x), so the boundary 
condition of fixed S Q appears to be very close to the 
boundary condition S — of Ref. [l[ . 

3. Correlation effects 

Finally we comment on the statement of Ref. 
that the statistical weight of a single vortex penetrat- 
ing through the film edge should be P ~ L/w, instead 
of P ~ L/£ used in It is noteworthy that the mod- 
els of Refs. [H-0| only hold in the limit of exponentially 
low density of vortices, thus P should coincide with its 
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value in the thermodynamic limit. The assumption of 
P ~ L/w is therefore inconsistent with the thermody- 
namics of vortices in thin films 0, El used to obtain 
P ~ L/£ in Ref. [3]. Here P ~ L/£ is the ID analog of 
the statistical weight P = C(L/£) 2 of a single vortex in 
the film of area L 2 where C ~ 1 depends on the distri- 
bution of the order parameter in the vortex core 13 1 . 

The assumption P ^ L/w resulted from the interac- 
tion radius ~ w of two vortices in the middle of the strip. 
However, uncorrelated hopping of vortices described by 
Eq. ^ imply that they enter the film at random times 
and are separated by distances larger than w at any 
given moment. Taking vortex correlations into account 
requires solving coupled equations for the higher order 
correlation functions which cannot be described by Eq. 
©. The authors of Ref. 0] selected rare events when two 
vortices enter the film nearly simultaneously and ascribed 
their statistical weight P ~ L/w to all vortex jumps. 
However, repulsion of vortices suppresses their simulta- 
neous entering the film, forcing them to go one by one so 
that a vortex can enter at any of L/£ edge sites after the 
preceding vortex in the area ~ w has already crossed the 
film. Such uncorrelated jumps Q have a much higher 
probability than the correlated jumps assumed in Refs. 

In addition, the interaction radius of vortices at 
the film edge (x ~ £) is much smaller than w because 
currents of two vortices spaced by the distance s along 
the edge are nearly extinguished by their antivortex im- 
ages, resulting in the dipole interaction U(£, s) ~ e(£/s) 2 
which does not extend well beyond s > £. 

In conclusion, we show the importance of the vortex 
core energy, the realistic behavior of U(x) at the film 
edge, and the physical boundary conditions to Eq. © 
for the calculation of thermally-activated hopping of vor- 
tices across narrow films. Disregarding these issues in 



Refs. P, Q has produced large numerical errors and a 
significant discrepancy between the results of Refs. 0, Q 
and [3j. This can also be essential for the interpretation 
of experimental data on thin film photon detectors. 
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